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ABSTRACT 

The cumulative size-frequency distributions of impact craters on planetary bodies in the solar 
system appear to approximate a universal inverse square power-law for small crater radii. In 
this article, we show how this distribution can be understood easily in terms of geometrical 
statistics, using a de-Sitter geometry of the configuration space of circles on the Euclidean 
plane and on the unit sphere. The effect of crater overlap is also considered. 
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1 INTRODUCTION 

When observations by the spacecraft Dawn enabled the first de- 
tailed survey of impact craters on the astroid 4 Vesta recently, 
it was noted that the cumulative size-frequency distribution can 
be approximated by a power-law with (double logarithmic) slope 
— 1.9 ± 0.1 for small craters, close to the geometric saturation slope 
of —2 l lMarchi et al. I d2012h ). This approximate inverse square 
power-law has also been observed on other pla netary bodies in the 
solar system, for example on Mars (see, e.g., IWerner & T anaka 
d20 1 lh and references therein). Since size-frequency distributions 
of impact craters are a powerful tool in planet ology to date surfaces 
and re construct planetary evolution (see, e.g., Michael & Neukuml 
d2010h . and the review bv lHartmann & Neukum I ( l200ll) for Mars), 
it is important to understand their origin. The fi rst detai l ed sta tis- 
tical model of impact cratering was given by M arcus ] dl964l) in 
the context of lunar exploration. Assuming that craters are circular 
but may overlap, are Poisson-distributed over the lunar surface sub- 
ject to a time-dependent production function with probability den- 
sity p(r, t) for a crater radius to lie in the interval (r, r + dr), he 
showed in a subsequent work that for any p(r, t) oc r _7_1 , 7 > 2, 
the cumulative frequency of craters larger than radius r is indeed 
approximately p roportional to r~ 2 in a steady-state at late times 
<M arcus 1 11966)). Therefore, when this state is reached, the distri- 
bution of craters is the same for a class of underlying production 
functions. Moreover, if a population of craters of radius r reaches 
geometric saturation, then its frequency is exact ly proportiona l to 
r~ because of circle packing on the surface (cf. lGault I fl970h ). It 
may also be recalled that power-law d istributions are self-similar 
(for craters see, e.g., lTakavasu I (l99G), pp. 18f and 34ff), but this 
feature alone does not explain the fairly consistent power of —2 
itself. 

The universality of this property, then, seems to lend itself 
naturally to a geometrical explanation. In this article, we present 
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a new geometrical approach to study crater distribution functions, 
which is inspired by ap plications of meas u res to cosmological mod- 
els of inflation (e.g., iGibbons & Turok ] J200i)). While any size- 
frequency distribution of craters is, in the absence of other resurfac- 
ing, the result of the production function which depends on prop- 
erties of the impactors creating the craters, it is possible that the 
size-frequency distribution evolves to a steady-state at least par- 
tially independent of its creation history, as shown by Marcus' re- 
sult. In such a case, one might expect a priori a uniform distribu- 
tion over the configuration space of craters. In order to investigate 
this point, we shall start in this paper with the assumption of a uni- 
form distribution of craters over their configuration space and show 
that this implies in fact an inverse square power-law for the size- 
frequency distribution of small craters. Of course, by contraposi- 
tion, this means that any deviation from this relationship implies 
a non-uniformity of the crater distribution over their configuration 
space, which in turn may have physically interesting implications 
for the underlying class of production functions. 

The structure of this paper is as follows. In Sec. [2] we con- 
sider the 3-dimensional configuration spaces of circles on the Eu- 
clidean plane and on the unit sphere, and it turns out that these 
can be given a de-Sitter geometry. By interpreting the configura- 
tion space of craters as that of circles and computing its canonical 
volume measure, we show in Sec.[3]how an inverse square power- 
law for the cumulative size-frequency distribution of small craters 
emerges. This model also shows that the global crater distribution 
is no longer self-similar once crater size becomes comparable to the 
size of the planetary body. Initially, we shall assume that craters are 
sparse with negligible overlap in the sense that their area filling fac- 
tor F< 1, corresponding to a small fraction of geometric satura- 
tion which has filling factor tt/(2\/3) « 0.907 (cf. lGault I ( fl97bT) ). 
This appears to be a reasonable assumption, as the observed dis- 
tribution of craters is typicall y below 10% o f geometric saturation 
(e.g. for the Moon, see again Gault (1970) who compares it to a 
"Mare Exemplum" created in laboratory experiments). The effect 
of crater overlap on the distribution function is then considered as 
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well. The conclusions are in Sec. [4] where an interpretation of this 
result and further applications are suggested. 



2 DE-SITTER CIRCLE GEOMETRY 

We begin by showing that the geometry of the configuration spaces 
of circles on the Euclidean plane and on the unit sphere in Eu- 
clidean space can be regarded as de-Sitter spacetimes. The first 
case will be seen in Sec. 12.11 as special case of the configuration 
spaces of n-spheres in n + 1 -dimensional Euclidean spaces. The 
discussion of this more general case will also help to establish the 
second case in Sec. 12.21 Then, in both cases, one can easily derive 
the canonical volume measures for these configuration spaces. The 
connection with more general sphere geometry is outlined in the 
Appendix. 



2.1 On the plane 

Consider n-dimensional unoriented spheres S" in the n + 1- 
dimensional Euclidean space E n+1 . Any such sphere is uniquely 
defined by its centre at a G E n+1 and its radius r > 0. Since 
these are independent of each other, the configuration space of such 
spheres M n E +2 = {§" C E n+1 } is n + 2-dimensional. The spe- 
cial case n — 1 corresponds, of course, to circles on the Euclidean 
plane. We will now see how the configuration space M E +2 can 
be interpreted as a n + 2-dimensional de-Sitter spacetime, and this 
spacetime of constant curvature can be described in the usual way 
as a quadric hypersurface in n + 3-dimensional Minkowski space- 
time E 1,n+2 with metric r] ab = diag(-l, 1, 1). 

Firstly, one may regard the Euclidean space E' I+1 as the inter- 
section of the null cone through the origin in an n + 3-dimensional 
Minkowski spacetime E 1,7l+2 with a null hypersurface, as follows. 
Given an arbitrary x G E n+ , define 



x 2 + 1 x 2 

— n — > x > — 



i \ g E l,n + 2 



where the square of vectors in E n+1 is with respect to the corre- 
sponding Euclidean metric, and let 

N a = (1,0, ...,0,1) G E 1,n+2 . 

Then each x is identified with a point which is both on the null cone 
through the origin, 

ria b X a X b = 0, 

and on the hypersurface 

Vab X a N b = -1, 

which is null since ri ab N a N b = 0. 

Secondly, consider a spacelike vector Y a = 
(Y°,Y,Y" n+2 ) G E 1 '™" 1 " 2 , so that 7] ab Y a Y b > 0, which is 
perpendicular to X a , 



r)a b X a Y h = 0. 



(1) 



Writing Y ■ N = r) ab Y a N b = -Y° + Y n+2 for short, one can 
recast condition |Q3 as 



x + 



Y ■ N 



y ab Y a Y b 
(Y • AT) 2 



(2) 



So for any k ^ 0, a given kY a uniquely defines an n-sphere in 
E n+ with centre and radius 



Y lvabY a Y» 
a = — — — — , r — 



Y ■ N' 



(Y ■ N) 2 ' 



respectively. Because of the freedom to rescale Y a we may set, 
without loss of generality, 



Vab Y a Y b = 1, 



(3) 



and use 10 to express Y a in terms of the variables defining the 
corresponding sphere, 



Y° = 
Y = 

yn + 2 



1 / a 2 + 1 



2 V r 
a 

5 

r 

1 /a 2 -l 



(4) 
(5) 
(6) 



Hence, this configuration space of spheres in Euclidean space is 
seen to be an n + 2-dimensional de-Sitter spacetime dS n+2 , 

M n E +2 = {Y a G E 1 '™" 1 " 2 : ri ab Y a Y b = 1} = dS n+2 , 

as required. The metric on this de-Sitter spacetime is induced from 
the Minkowski metric of E 1 '™" 1 " 2 , 



ds 2 = -(dV ) 2 + dY 2 + (dY 
= ^(-dr 2 +da 2 ), 



n+2\2 



(V) 



using <|4]-{6j . With yfj, as the metric of the de-Sitter spacetime, one 
can rewrite 10 as 

ds 2 = g® b dy a dy b , y a = (r,a) G dS n+2 , 
to obtain 

det g^drda 1 . . . da n+1 



dV, 



„n + 2 



drda 1 . . . da n+1 



(8) 



as the volume measure of the de-Sitter and hence the configuration 
space. Finally, as a special case of Eq. HJl one finds 



dVg = -TT drda 1 da 2 , 



(9) 



the volume measure of the 3-dimensional configuration space M E 
of circles on the Euclidean plane. 



2.2 On the sphere 

In view of the application to the distribution of craters on a plane- 
tary surface, we shall proceed by extending the previous case to the 
configuration space of circles S 1 on the 2-dimensional unit sphere 
S 2 centered at the origin of E 3 . Each circle is fully characterized 
by its angular radius p G (0, 7r/2) and by its centre n on the unit 

2- sphere in Euclidean 3-space, so that n 2 = 1 and the two spher- 
ical coordinates of the centre can be computed in the usual way. 
The corresponding configuration space M% is, as before, clearly 

3- dimensional and can also be considered a de-Sitter spacetime. 

To see this easily, one can utilize a trick to discuss this problem 
in terms of the previous one. Each circle on our unit 2-sphere § 2 
may be regarded as the intersection of another 2-sphere S 2 of radius 
r centred at a G E 3 with § 2 . As noted in the previous section, 
the configuration space of 2-spheres in Euclidean 3-space is 
a 4-dimensional de-Sitter spacetime defined by the quadric l[3} in 
5-dimensional Minkowski spacetime E 1 ' 4 with coordinates (|4}j6]l. 
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Now, for a given circle of angular radius p, there is a family of 
spheres § 2 whose intersection with the unit 2-sphere produces this 
circle. Hence, their radius r depends on the choice of a which, by 
symmetry, must be parallel to n. As shown in Fig. Q] (left panel), 
we can choose a such that r = tan p and a = n sec p. Then from 
<EH6jl, one finds Y 4 = and 



y u = - cot P , 

Y = — ncscp. 



(10) 
(11) 



With this choice, we have coordinates Z a = (Y°,Y) G E 1 ' 3 on 
the 4-dimensional Minkowski subspace Y 4 = which determine 
uniquely the centre and angular radius of a circle on the unit 2- 
sphere and, since ri ab Z a Z b = — (Y ) 2 + Y 2 = 1, define a 3- 
dimensional de-Sitter spacetime as the corresponding configuration 
space, 



Ms = {Z a G E 1 ' 3 : Vab Z a Z b = 1} = dS 3 , 

as promised. The metric on this de-Sitter spacetime is again in- 
duced by the ambient Minkowski space E 1 ' 3 , 

ds 2 = -(dY ) 2 + dY 2 = esc 2 p(-dp 2 + dn 2 ), 

from d 10141 U . subject ton 2 = 1. One can write the line element on 
the unit sphere in terms of the usual spherical coordinates (6, <f>), 

dn 2 = d<9 2 + sin 2 0d0 2 , 

and thus express the line element of the configuration space as 

ds 2 = gabdz a dz b , z a = (p, 9, <j>) G dS A , 
with metric 

(12) 



S { 2 2 2 2 \ 

9ab = diag ( — esc p, esc p, esc p sin 6) 



Thinking of a sphere in E 3 of radius Ro with a circle of radius 
R measured on the spherical surface so that its angular radius is 
p = R/Ro, then this metric clearly reduces to the one in Eq. 
10 in the limit of small circles, p <C 1, as expected. The volume 
measure corresponding to metric d!2t becomes 



dV§ 



•detfl&dpdfcty 



= esc 3 p sin 9dpd9d(f> = esc 3 pdpdui, 



(13) 



for the 3-dimensional configuration space M s of circles on the unit 
sphere, where dej = sin 6d8d(j) is the solid angle element on the 
unit sphere. 



3 APPLICATION TO CRATER DISTRIBUTIONS 

Having determined the geometry of the configuration space M % of 
circles on the Euclidean plane and M | of circles on the unit sphere, 
we shall now show how this may be applied to the distribution of 
craters on planetary surfaces. The assumptions of our model shall 
be that, firstly, each crater can be represented by a circle and, sec- 
ondly, the distribution of circles is uniform on their configuration 
space as defined in Sec. [2] The latter assumption therefore serves 
as a null hypothesis about the structure of the configuration space 
of craters. However, the configuration spaces of circles allow ar- 
bitrary overlap whereas, naturally, this is not the case for craters 
erasing older and smaller ones within. In Sec. 13.11 we shall con- 
sider the case that craters are sparse so that their overlap may be 
ignored. This condition is relaxed in Sec. 13.21 where we attempt to 
indicate the effect of overlap on the distribution function of craters. 




01 (p -p) 



Figure 1. Left: A circle of angular radius p at n on the unit 2-sphere S 2 
(shown as solid circle) in Euclidean 3-space is represented as the intersec- 
tion of a 2-sphere § 2 of radius r = tan p at a = n sec p (shown as dashed 
circle) with the unit 2-sphere. Right: The centre of a crater of angular ra- 
dius p wholly overlapped by a crater of angular radius p' may lie within the 
shaded area of solid angle u)(p' — p). 



3.1 Sparse craters 

3.1.1 On the plane 

Suppose that an infinitesimal area d^4 in the Euclidean plane con- 
tains the centres of dn craters of radii in the interval [r, r + dr] . 
Then the corresponding distribution function Je may be written 
dn = fEdrdA so that, given our assumptions and the result from 
Sec. 12. II dn oc dVj and hence 

Cp 

= 7# ( 14 > 

from Eq. ©, where Ce is a constant. If this constant is sufficiently 
small, the assumption of sparse craters will hold. In order to de- 
termine this and estimate its value, let ng be the total number of 
craters with radius r m i n r ^ r max wholly within a disk of ra- 
dius R ^ r max . Then the area filling factor of these craters is 



itr /s(r)dr = kCe In ■ 



F E 



from d 14b . so that for sparse craters with Fe < 1 we must have 



C E < 7T In 



(15) 



As R — > oo and the condition H5i continues to hold, tie/ it R 2 
tends to the total frequency of craters with this radial range per unit 
area on the Euclidean plane, 



fsdr, 
whence, using again 
Ce = 2vsr 2 niVL I 1 



2 

mill 



(16) 



The cumulative frequency u(r) of craters with radius greater than 
r per unit area as a function of size can also be obtained from the 
distribution function d 14b . 



,(r) = 



t { '\j ' Ce ( 1 1 
f E (r dr =—-—-_ 
2 V r 2 r 2 



(17) 



Finally, note that the approximation r m i n <JC r max <JC R applied to 
Eq.s d!54417b . which is reasonable for a crater count, gives Ce ~ 
2z/Er min <C 1 and an inverse square power-law for the cumulative 
frequency, 



/ \ _ mm 

v{r) ~ve — . 



(18) 
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3.1.2 On the sphere 

When the maximum radius of the craters is not negligible compared 
to the radius of the planetary body then, given our assumptions, the 
configuration space of circles on the unit sphere from Sec. 12.21 is 
applicable. Let Ro be the radius of the planetary body and R be 
the radius of a crater measured on the spherical surface so that its 
angular radius is p = R/Ro £ (0, 7r/2) as before. Now suppose 
dn is the number of the centres of craters with angular radius in the 
interval [p, p + dp] within the solid angle element duj = dA/Rg, 
so dn = fsdpduj with distribution function fs oc dVg given by 

fs(p) = C S esc 3 p, (19) 

from J 13t . where Cs is a constant. A crater of angular radius p 
subtends a solid angle 



(20) 



p2iv pp 

u (p) = / / srn dd8d(f> — 27r(l — cos p) 
Jo Jo 

on the planetary surface, so the area filling factor of craters with 
angular radius p m in ^ r ^ n/2 is then 



F s 



u(p)fs(p)dp 



, ■ i COS Omin Pmin 

nCs ( — h In cot — — 



1 + COS p m i n 

For the assumption of sparse craters to hold, we require that Fs <C 
1 as before. Now if p m in <C 1, that is, the smallest craters counted 
are much smaller than the radius of the planetary body, one must 
therefore require that 

Cs « (?t+k In— ) . (21) 

V 2 pmin / 

The total number of craters with radius greater than p mm is 

f 47T fi 

ns = fs{p)dpdu) 



o J p 



2nCs ( COt pmin CSC Pmin + In COt 



Pmin 



(22) 



Again if p m in <C 1, then we can approximate l !22t by the leading 
monomial to find 



C s 



n-SPmin 

2tv ' 



(23) 



The total cumulative number of craters with radius greater than p 
is also obtained from f 19b . 



n(p) 



Atv I fs(p'W 
' p 



= 2-kCs ( cot pese p + In cot ■ 



(24) 



This function is shown in Fig. [2] for three values of Cs, for which 
condition J2U holds in the range shown. Also, if p m m <§C 1, one 
can use the approximation l |23l > for the constant in Eq. J24b to find 

n(p) w nsp min (cotpesep + In cot . (25) 

This expression clearly recovers Eq. J 1 8 b in the limit of p <C 1 
where 4ttRqU = n and pRo = R ~ r, as expected. Likewise, one 
can easily convert d24t > to the cumulative crater frequency per unit 
area of the planetary body as a function of crater size. Note also 
that the global cumulative number J24t is given by trigonometric 
functions and not a power-law. Hence, the crater distribution is no 
longer self-similar, as is the case for small craters or for craters on 
the Euclidean plane in the limit r max — > co of Eq. d!7t . 




Figure 2. Double logarithmic plot of the total cumulative number n of 
craters larger than angular radius p under the assumption of sparse craters 
(solid curves), for p m in = 10 -3 and three values of the constant in Eq. 
{24): from top to bottom Cs = 10 -2 ,lCr 3 ,lCr 4 , respectively. Note that 
the curves approach the inverse square power-law n oc p~ 2 for small p. 
The dotted curve shows the effect of crater overlap for the top curve with 
A = —0.1 inEq. {32). 



3.2 Overlapping craters 

The spherical model discussed in Sec. l2.2l is. of course, the gener- 
alization of the planar model of Sec. 12. 1 l applicable to sparse craters 
on a planetary surface. However, a more realistic model for the 
distribution of craters should take into account the possibility of 
overlap. In this section, we shall consider how the more general 
spherical model can thus be extended. For simplicity we shall as- 
sume that, firstly, there is a finite period of crater formation and, 
secondly, the crater distribution on the planetary surface is always 
isotropic. Thirdly, we take l |19t to be the intrinsic distribution of 
craters produced, which may differ from the distribution observed 
on the planetary surface due to overlap. Finally, we assume that the 
centre and radius of a crater can be identified as long as a sector 
of its perimeter survives, that is, only craters that are wholly over- 
lapped by a later and larger crater are considered to be obliterated. 

Hence, one may take the initial distribution to be fs(p) and 
the final observed distribution fo(p) to be a function of p only. 
Craters of radius p will have been wholly overlapped by larger 
craters of radius p > p if their centres are within a solid angle 
of 



>{P ~ p) = 2tt(1 - cos(p' - p)) 



(26) 



from Eq. J20t . concentric with a larger crater of radius p' in the ob- 
served distribution, of which there are 4nfo{p')dp'. This is illus- 
trated in Fig.Q](right panel). In addition, there remain 47r/o (p)dp 
craters of radius p in the distribution observed on the planetary sur- 
face. By assumption, we take the distribution of craters of radius p 
that are produced to be proportional to fs(p)- Then this applies to 
the overlapped craters as well as to the sum of the observed craters 
plus the overlapped ones, of course with different constants of pro- 
portionality. We therefore expect to have 

fo(p) + f 2 fo(p')fs(pMp' - p)dp' oc f s (p). (27) 
J P 

It may be useful to express the distribution function as 

fo(p) = Cof( P ) 

where Co is the scaling factor proportional to the total number of 



© 0000 RAS, MNRAS 000.[T]{6] 



On de-Sitter Geometry in Crater Statistics 5 



observed craters no greater than p n 
limit of sparse craters, 



. which reduces to Cs in the 



(28) 



(29) 



no = 47rCo / /(p)dp. 

J Pmin 

Hence, using OH and l|26]l, one can rewrite as 

/(p) = !?(p) + A f 2 f(p')K(p,p')dp', 
where 

g(p) = esc 3 p, 
and 

P) = CSC 3 p(l - COs(p' - p)), 

and A < is a (by definition negative) constant. The problem of 
overlapping craters, at least in this simplified form, can therefore be 
expressed as the integral equation l |29l >, which is the standard form 
of a linear, inhomogeneous Volterra equation of the second kind. 
Note that its kernel is of separable form since it may be written 

K(fi,p) =^2ki(p)k'i(p), 

i 

where the ki and k' t are functions only of p and p', respectively. Al- 
ternatively, one can turn the integral equation l |29t into a differential 
equation. Letting 

Kp) = ^ 

we get 

- d f d 3 f 
dp dp 6 

a linear, homogeneous ordinary differential equation of third order 
with non-constant coefficients. Furthermore, one can define 



/GO 



din/ 



dp 

to reduce the order of ( 



i and convert it to 



Aff + f + f 



J dp + dp 2 



0, 



(3D 



a non-linear, inhomogeneous ordinary differential equation of sec- 
ond order with constant coefficients. While a detailed discussion of 
the properties of Eq.s J29M3U and the numerical comparison with 
observed crater distributions is beyond the scope of the present ar- 
ticle, we shall use the Volterra equation \29\ to illustrate the effect 
of little overlap for craters with small lower limit of the angular ra- 
dius, p m in <C 1. By this we mean a small deviation of fo from the 
sparse crater distribution so that the constant in the Volterra equa- 
tion is small, jAj <C 1. Then an approximate solution of Eq. j29\ 
for / can be expressed naturally as the leading terms of a Neumann 
series expansion of the full solution, 

oc 

/(p) = X> j ^o>), 

3=0 

whence the first approximation will be given by 

/(p) « u + Aui 

= esc 3 p + A / esc 3 p esc 3 p(l — cos(p' — p))dp' 
J P 

= csc 3 p(l + ^ f-cosp + lncot^) ] . (32) 



Integrating {32} in Eq. ( 128b . the leading monomial in p n 
the approximation 



yields 



no 



2iyCq 

Pmin 



4 V Pmin 



Hence, the constant A indicates the importance of crater overlap. 
In the current formulation, it is a free parameter that will depend 
on the total number of craters produced and on the number no of 
craters that remain observable. 



4 CONCLUSION 

Starting from the assumption that impact craters are uniformly dis- 
tributed over their de-Sitter configuration space, it is shown by Eq.s 
dl8t and d25b that this can indeed explain the approximate inverse 
square power-law for the cumulative size-frequency distribution of 
small craters, which is widely observed as noted in the Introduction. 
Moreover, while Eq. i24l shows that the global crater distribution is 
not self-similar, Fig.[2]illustrates that the deviation from the inverse 
square power-law due to the finite size of the planetary body is in 
fact mostly negligible, so that the crater distribution is effectively 
self-similar. Also, since smaller craters are relatively more affected 
by overlap, we expect that they are systematically undercounted 
relative to the distribution of sparse craters with negligible overlap. 
This is the case, as indicated by Eq. d32| l and the dotted curve in 
Fig.|2]for small p. 

Hence, the approach promoted in this article suggests that the 
de-Sitter configuration space may be a useful prior for the study 
of cratering histories. Deviations from the inverse square power- 
law correspond to non-uniform distributions of small craters on this 
configuration space. We suggest that such deviations from unifor- 
mity will indicate physically interesting processes like resurfacing 
due to igneous activity and erosion, as well as properties of the 
impactors before a steady-state of the crater population, as in the 
formalism of Marcus mentioned in the Introduction, is reached. A 
first step in this direction would be to classify production functions 
and resurfacing processes in terms of the deviation of the resulting 
crater configuration space from this de-Sitter geometry, which will 
depend on time. Hence, comparing these templates with observed 
maps of crater configuration spaces may yield new insights into 
cratering histories. 

Finally, this article notes a perhaps surprising formal connec- 
tion between the geometry of a cosmological spacetime and a prob- 
lem in planetology. Recall that the 3-dimensional configuration 
space of circles, and hence of craters as discussed, was interpreted 
as the 3-dimensional de-Sitter spacetime with the metric given by 
l l 1 2t and hence the line element, 



d2 2 i 2 , 2 

s = — esc pdp + esc 



p(d8 



2 +sin 2 6d(/) 2 ). 



Now in order to connect this with the de-Sitter spacetime familiar 
from cosmology, one can use the coordinate transformation esc p — 
cosh t so that 



ds 2 = — dt 2 + cosh 2 t(d9 2 + sin 2 



where t takes its usual interpretation as coordinate time. Extending 
our argument from two to three spatial dimensions, one obtains a 
measure on the space of spheres in 3-dimensional Euclidean space 
from Eq. l[8j for n — 2. This may be useful to study the statistics 
of bubbles in the universe, for instance supernova remnants within 
the interstellar medium or voids in the large-scale structure. Also, 
eternal inflation treats of bubble universes which are expected to 
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collide (see, e.g.. iGarriga. Gufh~& Vilenkin] d2007l) ). so finding a 
version of the overlap formula l |29t applicable to this problem may 
further elucidate the collision process of bubble universes. 
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APPENDIX 

Here we give a brief discussion of the geometry underlying the 
construction in Sec. [2] Circle geometry, beginning with the classi- 
cal problem of Apollonius to find the circles touching three given 
ones, is a special case of the sphere geometries developed by Lie 
(1872) and L aguerre (1881). Following partially the treatment by 
iBenz I d2012h . consider the Euclidean space E n+1 with its stan- 
dard Euclidean metric. Then a Laguerre cycle is either a point or 
an oriented n-sphere of radius r / centered at a £ E n+1 , with 
orientation expressed by the sign of r. Laguerre cycles can thus 
be regarded as elements of a real vector space V = R © E n+ . 
Assigning coordinates c = (r, a) to a Laguerre cycle is called cy- 
clographic projection, and addition of Laguerre cycle coordinates 
is defined as for vectors. The zero Laguerre cycle is = (0, 0). 
Now given two Laguerre cycles ci = (n, ai), C2 = (^2, a2), the 
real number 



P(ci,c 2 ) 



-(r 2 -ri) + (a 2 - ai) 



(33) 



where the product of vectors in E n+1 is understood to be the inner 
product with respect to the Euclidean metric, is called power of 
ci and C2. It turns out that ci and ci t ouch each other respecting 
orientation if, and only if, teenzl ( l2012h . proposition 3.43) 



and this is exactly the reflexive and symmetric contact relation stud- 
ied in sphere geometry. Moreover, V can be equipped with a prod- 
uct defined as 

C1C2 = — r\T2 + aia2. 

Letting dc = (dr, da), we therefore find that 



Ac 1 



-dr + da = P(c, c + dc) 



(34) 



measures the infinitesimal deviation of Laguerre cycles with re- 
spect to the contact relation, and this is seen to define a Minkowski 
metric on V . Note that J34b may be regarded as a line element on 
V which is invariant under transformations of Laguerre cycles pre- 
serving power J33l >. Hence, there is a close connection between 
the Minkowski spacetime of special relativity and the geometry 
of Laguerre cycles, w hich appears to have been pointed out first 
by iTimerdingl jl912h . In addition to the set of Laguerre cycles T, 
one can introduce the set of oriented Euclidean hyperplanes E in 
E n+1 and the object infinity 00, and define corresponding con- 
tact relations. Then A = r U E U {00} is called the set of Lie 
cycles, and bijections A : A — > A that preserve contact rela- 
tions form t he group of Lie transformations of E n+1 . It turns out 
jBenz I ^20121) . proposition 3.56) that to every Lie cycle a homo- 
geneous Lie cycle coordinate L — [L° , L, L n+2 , L n+Z \, where 
[L°,L,L n+2 ,L n+3 ] = [kL°,kL,kL n+2 ,kL n+3 ] for any k ^ 0, 
can be assigned bijectively, which satisfies the Lie quadric 

-(L°) 2 +L 2 + (L" +2 ) 2 -(L" +3 ) 2 



0. 



(35) 



For a Lague rre cycle c = (r, a), the Lie cycle coordinate is (c.f. 
IBenz 1 {2012), p. 154, applying a sign change and reordering) 



P(0,c) + 1 P(0,c)-1 



L(c) = 



By homogeneity of the Lie cycle coordinates, we can set L n+3 
to obtain 



L(c) = 



a-' + l 



r 

n+2 



-r , — ,-- 



using J33b and the coordinates {4]l-l[6]l defined in Sec. [2] Then these 
are seen to satisfy 

j-yn + 2^2 



(y°) 2 + Y 2 



1 



(36) 



because of the Lie quadric ( 135b . Suppose now that we consider the 
line element 



ds 2 = 



-(dr 



0n2 +dY 2 + (dF n+2 ) 2 



(37) 



1 I A 2 , a 2\ dc 



using OU. Clearly, any transformation (Y° , Y, Y n+2 ) M> (Y° + 
dY°,Y + dY,y n+2 + dy n+2 ) that preserves our choice (gUl 
leaves this line element invariant, and in fact null ds 2 = so that 
dc 2 = 0. Hence, such transformations preserve the contact rela- 
tions of the corresponding Laguerre cycles, as described above. 
Under this assumption, then, (Y° ,Y ,Y n+2 ) may be regarded as 
a point on the de-Sitter configuration space of Sec.|2]defined by the 
quadric hypersurface d36b in the Minkowski spacetime E 1 > n + 2 with 
line element d37b - A connection between Laguer re sphere geo metry 
and de-Sitter spacetime seems to emerge first in iGraf I I I 1934b . who 
uses the opposite metric signature. 



P(ci,c 2 ) = 0, 
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